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Abstract 

We consider a generalisation of the classical Lehmer problem about 
the distribution of modular inverses in arithmetic progression, intro- 
duced by E. Alkan, F. Stan and A. Zaharescu. Using bounds of sums 
of multiplicative characters instead of traditionally applied to this kind 
of problem Kloosterman sums, we improve their results in several di- 
rections. 

1 Introduction 

Given modulus q > 2, we denote by U q the set 

U q = {n : 1 < n < q, gcd(n, q) = 1}. 

that is, j^lA q = (p(q), the Euler function. 

For n G U q we use n to denote the modular inverse of n, that is, nn = 1 
(mod q), n G U q . 

The classical question of D. H. Lehmer (see Problem F12]) about the 
joint distribution of the parity of n and n has been solved by W. Zhang 



Recently this question has been generalised by E. Alkan, F. Stan and 
A. Zaharescu P as follows. Given vector a = (oi, . . . , a^ + i ) G U^ +l and 
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b = (£>!, . . . , £>fc+i) G Z + we consider the set jV(a, b, q) of vectors n = 
(ni, . . . , n k ) G such that 

rii = bi (mod aj), i — l,...,k, 
ni...n k = b k+1 (mod a k +i). 

Generalising several previous results of various authors, (for instance, 
of 0H2H2DII2I1), E. Alkan, F. Stan and A. Zaharescu P have shown that 
for any fixed k, the bound 

#AT(a, b, q) = + O^- 1 '^) (1) 

a l ■ ■ ■ Q-k+l 

holds uniformly over all vectors a G U k+l and b G Z fc+1 . In particular, since 

S C loglog(, + 2) (2) 

for an absolute constant c > 0, see |T7J Section 1.5.4], we see that the 
bound is nontrivial for 

a x ...a k+1 <q x l 2 - 5 (3) 

for any fixed 8 > 0, provided that q is large enough. 

The main tool of [Q are bounds of exponential sum, in particular, mul- 
tidimensional Kloosterman sums. Here we show that using bounds of mul- 
tiplicative character sums, such the classical Polya- Vinogradov and Burgess 
bounds, see [TTH Theorems 12.5 and 12.6], one can improve the bound (JIJ. 
For example, we obtain a bound, which in particular implies that 

m(a) k /llall<7*~ 1+ °W a (k+i)/2+o(i) \ 

#A/-(a,b,g)= m) + 0[™- + wa+i)), W 

a 1 ...a k+1 V a i--- a fc+i (a! . . .a k+1 )W +1 ) J 

where ||a|| is the Euclidean norm of a, which is equivalent to (0) for k = 2 
and ||a|| = 0(1) and always improves it if either ||a|| grows together with q 
or if k > 3 (in this case, with respect to both dependence on q and a). 

We note that instead of (JHJ), the bound of our Theorem |H1 is nontrivial 
when simultaneously 

q (k 2 -i)/2k-s if 2 < A; < 4, 
|a|| < q 1 - 5 and a x ... a k+1 < <j q 5 / 2 ~ s if k = 5, (5) 

2(fc 2 -l)/3(fc+2)-<5 jf ^ > g 
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for any fixed 5 > 0, provided that q is large enough. In fact we consider 
a more general case when n%, . . . ,1%). and n\. . .n k belong to a certain box 
inside of the cube T^+i, where 

T s = (R/Z) s = [0, l) s 

is the fc-dimensional unit torus. 

The question about the distribution of elements of jV(a, b, q) in various 
regions of has also been studied in For an arbitrary region Q C we 
denote by N"a(a, b, q), the set of vectors n £ A/"(a, b, q) which belong to the 
dilated region qfl. Let A(J?) denote the Lebesgue measure of i?. It has been 
show in lj that for any fixed k and region i? C T^, with piecewise smooth 
boundary, 

#A^(a, b, q) = \(f2) ^ + 0( ? *-W+i)+ (i)) (6) 

a\ . . . dk+i 

holds uniformly over all vectors a £ U^ +1 and b £ Z fc+1 . 

Here we show that using some bounds from pQ in a combination with some 
results of M. Laczkovich ^3] and of H. Niederreiter and J. M. Wills [TS], from 
the theory of uniformly distributed sequences leads to a better error term in 
the asymptotic formula ©• 

Furthermore, we also consider a generalisation to the joint distribution 
of ni, . . . ,rik and rii . . . in arbitrary regions. Namely given an arbitrary 
region C T^+i we estimate the cardinality of A4e(sL, b, q), which is the set 
of vectors (n 1; . . . , n k ) £ A/"(a, b, q) for which (n 1; . . . , n k , n\ . . . n k ) belongs 
to the dilated region qO. 

Finally, in the case of prime q = p, we show that a result of A. Ayyad, 
T. Cochrane and Z. Zheng Theorem 2] leads to some improvements. 

We conclude with a short discussion of possible ways to improve our 
results and of some open problems. 

Throughout the paper, the implied constants in the symbols l O\ and 
'<C' may depend on integer parameters k and r and a region Q C T^. We 
recall that the notations U = 0(V) and V U are both equivalent to the 
assertion that the inequality \U\ < cV holds for some constant c > 0. 
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2 Preparations 
2.1 Character Sums 

Let X q be the set of all j^X q = (p(q) multiplicative characters of q. We refer 
to ^2] for definitions and basic properties of multiplicative characters such 
as x( u ) = f° r an y X £ if gcd(M, g) > 1 In particular, we recall that for 
u E Z, 



see [T51 Theorem 5.4]. We also use Xo to denote the principal character. 
In particular, we immediately see the following bound 

Lemma 1. For any integers K and L with 0<K<K + L<q, an integer 
a > 1 with gcd(a, q) = 1 and an arbitrary integer b, 



Proof. We recall that if gcd(n, q) = 1 for the conjugated character \ we have 
x(n) = x (n) ■ Therefore 
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£ E *(n) <¥>(?)(£/<*+!)• 



n=b (mod a) 



2 



2 









_ft:+i<n<ft:+L 

?i=6 (mod a) 



X^Xq K+l<n<K+L 
n=b (mod a) 
gcd(n,(/)=l 





n=fe (mod a) m=6 (mod a) 
gcd(n,q)=l gcd(m,,g)=l 






K+l<n<K+L K+l<m<K+L x£X q 
n=b (mod a) m=b (mod a) 
gcd(n,q)=l gcd(m,q)=l 
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Clearly n takes at most L/a+1 possible values and since < K < K + L < q, 
for each n the value of m is uniquely defined. Therefore T < L/a + 1, which 
concludes the proof. □ 

The following result is a combination of the Polya- Vinogradov bound (for 
r = 1) and Burgess (for r > 2) bounds, see [TUl Theorems 12.5 and 12.6]. 

Lemma 2. For any positive integers U and V < q, the bound 

u+v 

*m ^ v i - i/r v r+i)/4r2+o(i) . 

n=U+l 

holds with r = 1, 2, 3 for any q and with arbitrary integer r if q = p is prime. 
Now, using the identity 

_ftT+l<am+f><i : f+L 

(8) 



max 

X&Xq 

Xt^XO 



if+l<n<lf +L 
n=6 (mod a) 



(K+l-b)/a<m<(K+L-b)/a 

we derive from Lemma 121 the following useful estimate. 

Lemma 3. For any positive integers K , L and a > L such that gcd(a, q) = 1 
and an arbitrary integer b, the bound 



max 

Xt^XO 



K+l<n<K+L 
n=b (mod a) 



< g (4r 2 -3r+l)/4r 2 +o(l) a -(r-l)/r 



Zio/ds wt/i r = 1, 2, 3 /or an?/ q and with arbitrary integer r is q = p is prime. 

We also need to approximate the value of the sum of Lemma El with the 
principal character xo- 

We denote by w(q) the number of prime divisors of q. 

Lemma 4. For any positive integers K and L, an integer a > 1 such that 
gcd(a, q) — 1 and an arbitrary integer b, 

(p(q)L 



E 



K+l<n<K+L 
n=b (mod a) 



aq 



+ 0(2?®). 
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Proof. Clearly 



E 



E 



1. 



K+l<n<K+L 
n=b (mod a) 



K+l<am+b<K+L 
gcd(am+b,q)=l 



Using the Mobius function fi(d) over the divisors of q to detect the co- 
primality condition and interchanging the order of summation, we obtain 

E i = E^(| + o(i)) = §E^+o|Ewi 



if+l<am+&<,ff+,L 
gcd(am+&,(/)=l 



from which the result follows immediately. 



□ 



Finally, if q = p then we also use the following bound which follows from 
a result of A. Ayyad, T. Cochrane and Z. Zheng [21 Theorem 2] and the 
identity ©. 

Lemma 5. Let q = p be prime. For any integers K and L with < K < 
K + L < p, an integer a > 1 with gcd(a,p) = 1 and an arbitrary integer b, 



E 

X^XO 



K+l<n<K+L 
n=b (mod a) 



<p(L/a + l) 2 (logpf 



2.2 Discrepancy 

For a finite set T C T s of the unit s-dimensional set, we define its discrepancy 
with respect to a domain S C T s as 



-A(: 



where, as before, A is the Lebesgue measure on T s 
We now define the discrepancy of T as 

D{T) = sup A{F,LT), 

ncj s 



where the supremum is taken over all boxes 77 = [a±, /3i) x . . . x [a s , /3 S ). 
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As usual, we define the distance between a vector u G T s and a set r C T s 

by 

distfu, J 1 ) = inf llu — wll, 

where, as before, ||v|| denotes the Euclidean norm of v. Given e > and a 
domain S C T s we define the sets 

S+ = {u G T S \S | dist(u, S) < e} 

and 

S e " = {ueS| dist(u, T a \S) < e} . 

Let be an arbitrary increasing function defined for e > and such 
that 

lim hie) = 0. 

e^O 

As in [TH IT6] , we define the class 5^ of domains S C T s for which 

A (»+) < /i(e) and A (~7) < 

A relation between D(T) and ^(J? 7 , H") for S G 5^ is given by the follow- 
ing inequality of M. Laczkovich (see also [To]). 

Lemma 6. For any domain S G Sh, we have 

A^ 1 E)<^h{s 1 l 2 D{r) 1 l s ). 

Finally, the following bound, which is a partial case of a more general 
result of H. Weyl [IS] shows that if S has a piecewise smooth boundary that 
S G Sh for some linear function h(e) = Ce. 

Lemma 7. For any domain S G Sh with piecewise smooth boundary, we 
have 

\(~±)=0(e). 



3 Main Results 

3.1 Distribution in Boxes 

Here we study A4 e (a, b, q) in the case where = £ is a box E C T fc+1 and 
in particular we generalise and improve the bound (JTJ). 

We recall that we use ||a|| to denote the Euclidean norm of a G U\ +x ■ 
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Theorem 8. For r = 1,2, 3 ; any fixed k > 2, and a box 

E = [ax,Px) x ... x [a k+1 ,p k+1 ) C T k+1 

the bound 

ip(q) k 



#Ms(a., b, q) = X(E)- 



a\ . . . a k +i 

fe-l+o(l) fc-(3r-l)(/fc-l)/4r 2 +o(l) 



+o I ■ |a||r — + 9 



/io/ds uniformly over all vectors a G and b G Z fc+1 . 

Proof. We see from (JZj) that 

#Ms(a,b,g) = J2 ••• E 

aig<ni</3iij afc+i<?<nfc+i</5fc+i<? 
ni~h (mod ai) n fe+1 =6 fe+1 (mod a k+1 ) 

J-^^ x ( ni ...n fe+1 ). 

We now change the order of summation and note that by Lemma 0] the 
term corresponding to the principal character x — Xo is equal to 



E - E 



09(0) 

aig<ni</3iij a fe+ ig<n fc+ i</3 fe+ i<j 

m=bi (mod ai) n fc+1 =6 fc+1 (mod a k+1 ) 

gcd(ni,g)=l gcd(n fe+1 ,g)=l 



iff((M* + (^)) 

A(E) ^ + O f l*^2M^) 
ai . . . afc + i \ ai . . . a^+i 

A(r) y(g) fc ] c fMq k - 1+o{1) 



since 



a\ . . . a k+ i \ a*i . . . a k+ i 

i x „ log? 
log log q 
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see [T71 Section 1.5.3], and the bound (J2J). Hence, 



#M g a, b, g) = A( ^ +O lM ' 



Q>1 ■ ■ ■ Ctk+1 



a\ . . . afc + i 



+ R 



(9) 



where 
R 



— V T ••• V 
09(0) ^— ' ^— ' ^— ' 

XG^9 aiq<ni</3ig a fe+ i<?<n fe+ i</3j. +1 g 
XT^XO "1=1)1 (mod 01) n fe+ i=ife fe+ i (mod a fe+1 ) 



X (ni . . . Tifc+i) 



fe+i 



tpfq) YTl Y 

X^X q u=l a u q<n v <(3 u q 
X^XO n v =b v (mod a„) 



x M ■ 

Thus using the Holder inequality we obtain 



R k+1 < 



^ fe+1 

nE 



X+XO 



Y * 



a v q<n v <f3 v q 
n v =b v (mod a„) 



fe+1 



(10) 



Since the bound is trivial for ||a| > q, we can assume that 

max{ai, . . . ,a k+ i} < q. 

Applying Lemma HU to the (k — l)th power of the character sums for 
each u = 1, . . . , k + 1, and then extending the summation over all characters 
X G X, we obtain that for r = 1, 2, 3 



E 

X^Xq 

x¥=xo 



Y 



ot v q<n v <fl v q 
n v =b v (mod a„) 



< I /Y (4r 2 -3r+l)/4r 2 +o(l) /1 -(r-l)/r 



fc-1 



E 

xex q 



E 



q<n v <P 
n u =b u (mod a u ) 
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We now use Lemma d which implies 

fc+i 



E 

xex q 



E 



X\n V/ 



a v q<n v <f} v q 
n„=b„ (mod a„) 



< ( g (4r 2 -3r+l)/4r 2 +o(l) a -(r-l)/r 



fc-1 



q 2 a v x 



<T „(fe+l)-(3r-l)(fc-l)/4r 2 +o(l) -(rfc-fc+l)/r 



Substituting this bound in (fTUj) and using (j2J), we obtain 

-(rk-k+l)/r(k+l) 

^ < ?M 3r-l)(*-l)/4r 2 +0 (l) ( JJa/ 



>,i/=l 



which together with © completes the proof. 



□ 



In particular, taking r = 3 we see that the bound of Theorem |H1 implies 
that for any fixed k and 5 > there exists rj > such that under the 
conditions (0) we have 

Wa, b, g) = (A(E) + 0(g-")) 

ai . . . afe + i 

Moreover, using the trivial bounds 

||a|| < oi . . .a fc+ i and (ai . . . a fc+ i) 1/(fc+1) > 1, 
we derive from Theorem |H] that 



#M s ( a ,h,q) = \(E)- 



+ O (g fe - 1 +°( 1 ) + q (k+l)/2+o(l)^ _ 



ai . . . afc +1 

Finally, taking E = T^ +1 and r = 1 in Theorem |H1 we obtain (@J). 
3.2 Distribution in General Regions 

Here we give an improvement and generalisation of the asymptotic for- 
mula (El)- 
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Theorem 9. For r = 1,2,3, any fixed k > 2 and region O C T^ +1 with 
piecewise smooth boundary, 

#.M (a,b,g) = A(0)- 



O-l ■ ■ ■ a k+l 

+ O 



a || l/(fc+l)gfc-l/(fc+l)+o(l) ^fc-(3r-l)(fc-l)/4r 2 (fc+l)+o(l) 



ai...fl fc +i (ai.^k+i) 1 -^- 1 ^ 

ao/ds uniformly over all vectors a G and b G Z fc+1 . 

Proof. It follows from Theorem |H] and the estimates (J2J) and (@J) that we have 
the bound 

£>(.A(a, b, q)) « HaUg- 1 ^ 1 ) + ( fll . . . afc+1 )(r+*-l)/r(*+l) g -(3r-l)(*-l)/4r»+o(l) 

on the box discrepancy of the set 



.A(a, b, g) = <M — — , Ul ' ' ' Uk ) : (m, . . . , n fc ) G A/"(a, b, g) 

Therefore, by Lemma El and Lemma [7| we conclude that the discrepancy 
of A(a, b, g) with respect to O satisfies 

A{A{&,b,q),G) « || a ||i/(*+i)g-i/(fc+i)+o(i) 

+ (Oi . . . aA , +1 ) fe A(fc+l) 2 g -(3r-l)(fc-l)/4r 2 (fc+l)+o(l) 

which is equivalent to the desired result. □ 

Certainly applying Theorem El with — Q x [0, 1) where Q C one 
immediately obtains an asymptotic formula for 7Vj?(a, b, g) (which is already 
stronger than However since the problem of estimating A/^(a, b, g) is of 
lower dimension (A; instead of k + 1) one obtains a slightly stronger bound in 
this case. 

Theorem 10. For r = 1,2, 3, any /zxed > and region CTj mi/i piece- 
wise smooth boundary, 



#A/^(a,b,g) = X(f2)- 



ai . . . afc+i 
+0 



„ II l/k fc-l/fc+o(l) fc-(3r-l)(fe-l)/4r 2 fc+o(l) 

II y y 



a!...a fe+1 (o 1 ...a w . 1 )M^-DM(*+i) 
aoZds uniformly over all vectors a G and b G . 
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Proof. Taking the box E = II x [0, 1), where a box 77 C we see that 
it follows from Theorem |S] and the estimates (J2J) and (@J) that we have the 
bound 



D(£(a,b,g)) < ||a||g- 1+0 W + (en . . . ak+1 )(r+k-i)Kk+i) q -(3r-i)(k-iy4r-+o(i) 
on the box discrepancy of the set 

B(a, b, g) = | . . . , y ^ : (m,...,n k ) e jV(a, b, g) 

that is, of the same strength as that for the set A(a, b, q) in the proof of 
Theorem |UJ 

Therefore, by Lemma Inland Lemma [7| we conclude that the discrepancy 
of £>(a, b, q) with respect to Q satisfies 



< ||a|| 1/,6 g _x/ ^" w + (ai . . . afc+i^' ^"^ ' ^"^^" 
which is equivalent to the desired result. □ 



|l/fc -l/fe+o(l) , / \(r+fc-l)/rfc(fc+l) -(3r-l)(fe-l)/4r 2 fc+o(l) 



We remark that although in the case of = 2 and fixed a, Theorem |H] 
(with the optimal choice of r = 1) is equivalent to the bound of Theo- 
rem^] still improves © due to our use of Lemma El instead of the arguments 
from PP. 

3.3 Some Improvements for Prime q = p 

Here we show that if q = p is prime and k > 3 then using Lemma El instead 
of Lemma [U leads to a stronger bounds wit respect to the product a\ . . . cik+i. 

Theorem 11. Let q = p be prime. For any fixed integer r > 1 and k > 3, 
and a box 

E = [«!, ft) x . . . x [ajfe+i, C T fc+1 

£/ie bound 

#M E {*,b,p) = \{E)- '' 



+ O 



a\ . . . a^+i 

| a ||pfc-l+o(l) pfc-(3r-l)(fc-3)/4r 2 +o(l) 



ai...a fc+1 (a 1 ...a fe+1 ) 1 -( fe + 2 ''-3)/Kfe+i) 
fro/ds uniformly over all vectors a e a^d b G Z fc+1 . 
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Proof. We proceed as in the proof of Theorem |S] and also note that one can 
change (p — l) fc to p k in © without changing the error term. In particular, 
we still have (jlOj) . Again as in the proof of Theorem |H] we apply Lemma El 
however this time to the (k — 3)th power of the character sums for each 
v = 1, . . . , k + 1 (and this time we do not extend the summation over all 
characters x £ <^)> we obtain that for any integer r > 1 



E 

Xv^XO 



Oil/ 

n v =b v (mod a„) 



fc+1 



< 



P 



(4r 2 -3r+l)/4r 2 +o(l) _(r-l)/r 



fc-3 



E 

Xt^XO 



E 



X\n h 



a v p<n v <fi v p 
n v =b v (mod a„) 



We now use Lemma which implies 

fc+i 



E 

xg*p 

Xt^XO 



E 



a v p<n v <f3 l/ p 
n v =b v (mod a v ) 



< ( J 4?,2 -3>"+l)/4r 2 +o(l) -(r-l)/r 



fc-3 



.,4-2 



< ^+l~(3r-l)(fc-3)/4r 2 +o(l) -(rfc-fc-r+3)/r(fc+l) 

— v a v 



Substituting this bound in (fTU|) and using (J2J), we obtain 

/fc+i 

R < p fc -(3»-3)(fc-l)/4r 2 +o(l) J j-J 



fc+1 \ -(rfc-fc-f+3)/r(fc+l) 
u=l 



which together with (jUJ) completes the proof. 



□ 



In particular we see that the bound of Theorem ^2 taken with r = 1 
implies that for any fixed k and 5 > there exists rj > such that under the 
conditions 



\a\\<p l - s and a x . . . a fc+1 < j9 (fe+1)/2 - 5 



we have 



#M s ( a ,h,p)=(\(E)+0(p-v)) 



V 
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However, taking a sufficiently large r we see from Theorem ^2 that for any 
5 > there exists K and rj > such that for k > K the above bound holds 
under the condition 

|| a|| < p 1 - 5 and a x . . . a k+l < p^l^ k 

We also have analogues of Theorems 01 and [TU1 

Theorem 12. Let q = p be prime. For any fixed integer r > 1 and > 3, 
and region O C T^+i wt/i piecewise smooth boundary, 

#M e (a,b,j ) ) = A(0)- /r 



ai . . . afc+i 

I „ II l/(fc+l) fc-l/(fe+l)+o(l) fc-(3r-l)(fe-3)/4r 2 (fc+l)+o(l) 

+ 01^ ^ + 1 



ai ...a k+1 ( ai ...a k+i y-^-^i^ 2 

holds uniformly over all vectors a G U k+1 and b G Z fc+1 . 

Theorem 13. Let q = p be prime. For any fixed integer r > 1 and A; > 3, 

and region Q C «n£/i piecewise smooth boundary, 

#Afo(a,b,p)=A(tf)- /r 



ai . . . afc + i 

I „ II l/fc fc-l/fc+o(l) fc-(3r-l)(fc-3)/4r 2 A;+o(l) 

+ O I + ; 



ai...a fc+1 ( ai ...a fc+1 )MM*-3)/rk(fc+^ 
/io/ds uniformly over all vectors a G ZY^ +1 and b G Z fc+1 . 



4 Concluding Remarks 
4.1 Further Improvements 

Clearly, if some of a 1; . . . , a k +i are of different order magnitude, then in the 
proofs of Theorems |H1 and ^2 one can use Lemma [21 with various values of 
r for each v which may lead to stronger bounds. However it seems that 
the optimal strategy of applying these results heavily depends on various 
relations between the sizes of a\, . . . , a k +i and q. 
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We believe that there are several further possibilities of improving Theo- 
rem |H1 and in particular improving the threshold ©. Certainly there should 
be a variant of the result of A. Ayyad, T. Cochrane and Z. Zheng [21 Theo- 
rem 2], given in Lemma which holds for arbitrary composite moduli q (see 
also [4j). In fact, J. B. Friedlander and H. Iwaniec [7j give such a bound, but 
only for special intervals (starting at the origin). Certainly, obtaining such 
a general result is of independent interest. Furthermore, the technique used 
by M. Z. Garaev jH] can probably be useful as well. 

We also note that by a result of H. Niederreiter and J. M. Wills ^1 for 
the class of convex sets and Q the implied constants do not depend on the 
set (see also [5, Theorem 1.12] and Theorem 1.6, Chapter 2]). 



4.2 Some Open Problems 

It is certainly interesting to study various geometric properties of the set 
A/"(a, b, q). For example, let 

if (a, b, q) = max min |rij — n{ . . . rik\ ■ 

(ni,...n k )£j\f(a,h,q) l<i<k 

For k = 1, a = (1, 1), b = (0, 0), the value of 

H(q) = max \n — n\ 

n<=U q 

has been studied in [HI EJ- in particular, it has been shown in ^T] that 
H{q) = q + (g 3 / 4+ °«). It has also been shown in |Hj that H(q) is influenced 
by the distribution of divisors of qs — 1 for small values of s, and thus some 
lower bounds on H(q) have been derived. It would be interesting to find 
out whether the behavior of if (a, b, q) is also influenced by some arithmetic 
properties of the modulus q. 

Finally, one can also study various geometric properties of the convex 
closure of A/"(a, b, q). For example, for k — 1, a = (1, 1), b = (0, 0), that is, 
for the set 

N{q) = {{n,n) : n E U q } 

some lower and upper bounds on the number of vertices V(q) of its convex 
closure have been given in ^2]- These bounds as well as some numerical 
calculations suggest that the convex closure of Af(q) does not behave as 
the convex closure of a random set, but rather is affected by the arithmetic 



15 



structure of q — 1 (and probably of qs — 1 for small integers s). It would be 
interesting to see whether the same effect appears in the behaviour of the 
convex closure of jV(a, b, q) for larger values of k and "generic" vectors a 
and b. 
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